NAVAL  POSTGRADUATE  SCHOOL 
Monterey,  California 


© 


DTIC 


_ THESIS _ 

CROSS-SHORE  SEDIMENT  TRANSPORT  ON  A 
NATURALLY  BARRED  BEACH 

by 

Randall  T.  Humiston 
March  1993 

Thesis  Advisor  Edward  B.  Thornton 

Second  Reader  Thomas  C.  Lippmann 

Approved  for  public  release;  distribution  is  unlimited. 


>3  5 


93-12022 


m  nil 

i 


REPORT  DOCUMENTATION  PAGE 


form  Approved 
OMB  No  Q7C4-Q188 


0'vDl't  'eoorr.nq  Owr0«n  'af  I'M  iO'W'On  of  rt»orrt-»jt,on  -S  tO  4**'«q*  ‘  '‘■3ur  Off  fspot'Sf  «<  *<1**9  Iff  »-*R*  '#*•**»  *  *9  Mtr *,rt;C**  vfirt*  '•q  *»  s*  ^9  34.4 

9*thff'nq  *n<J  m4.nt41r.rtq  ?**  34t4  4rtd  ca'"0t*tirtt)  *r>a  ffvfA.rq  Thf  ■  :lifCTiCrt  3»  .rt?a','"4TtO*  Se«0  ECmrpfrtfi  ffOifflirq  tb-S  Ov»f 3f rt  Vm  ^*tf  I'  if-,  ;r».ff  4vOf<t  0»  **'  v 

:ci*fd»on  of  »nf0fm4t»o*v  «c«wj*ng  jvgg**t>on»  ^or  'Moong  tn.s  ourof"  to  Afasb-rtgiOH  xfaaaua'ters  Sfyces  C  rec*o'#!r  »of  Ooe'*t  c*s  j"<3  ns  ,  >5  **Vsv" 

0*vS  Krqhwiy.  Sw*t*  UO*  Arlington.  ./A  22202^302  4rtd  to  t*f  Q^-ff  O*  4n<3  9v,0qfT  »iOffwo^»  •fdgrt.O^  P'O.fd  ?0.’C*  0’t8>  A4i\rt"rtqt0rt  :  C  ;CSCJ 


1.  AGENCY  USE  ONLY  (Lave  blank) 


3  REPORT  TYPE  AND  OATES  COVERED 


4.  TITLE  AND  SUBTITLE 

CROSS-SHORE  SEDIMENT  TRANSPORT  ON  A  NATURALLY 
BARRED  BEACH 


6.  AUTHOR(S) 

Humiston,  Randall  T. 


7  PERFORMING  ORGANIZATION  NAME(S)  AND  AOORESS(ES) 

Naval  Postgraduate  School 
Monterey,  CA  93943-5000 


8  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


9.  SPONSORING  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


10  SPONSORING  MONITORING 
AGENCY  REPORT  NUMBER 


tl.  SUPPLEMENTARY  NOTES 


The  views  expressed  in  this  thesis  are  those  of  the  author  and  do  not  reflect  the  official  policy  or  position  of  the 
Department  of  Defense  or  the  U.S.  Government. 


12).  DISTRIBUTION  /  AVAILABILITY  STATEMENT 


12b  DISTRIBUTION  CODE 


Approved  for  public  release,  distribution  is  unlimited 


13.  ABSTRACT  (Maximum  200  worch) 


Bailard’s  sediment  transport  model  (1981)  is  evaluated  using  field  data  obtained  on  a  naturally 
barred  beach.  Principal  field  measurements  consist  of  a  cross-shore  array  of  bi-directional  current 
meters  spanning  the  surf  zone  and  daily  bathymetric  surveys.  The  model  predicts  bed  and  suspended 
load  transport  separately  based  on  various  velocity  moments.  The  velocities  are  partitioned  into  mean 
currents,  long  waves  (<  0.05  Hz)  and  short  waves  (>  0.05  Hz)  to  determine  their  relative  importance 
to  the  transport.  Velocity  moments  are  then  computed  over  90  minute  intervals  to  resolve  tidal 
fluctuations.  Finally,  predicted  transport  rates  are  integrated  and  compared  with  daily  cross-shore 
bathymetric  profiles  (averaged  over  a  400m  length  of  beach). 

Results  indicate  that  suspended  load  was  consistently  greater  than  bed  load,  as  much  as  an 
order  of  magnitude  during  episodes  of  large  incident  waves,  owing  to  the  slow  fall  velocity  (2cm/s)  of 
the  fine  grain  sand  within  the  surf  zone.  The  contribution  by  the  mean  current,  long  and  short  waves 
to  the  cross-shore  transport  were  of  the  same  order.  Variance  of  transport  during  all  stages  of  tide  and 
over  a  range  of  incident  waves  were  consistently  greater  in  the  vicinity  of  the  bar  and  trough  than 


1  A.  SUBJECT  TERMS 


Sediment  transport,  DELILAH,  Cross-shore  transport 


15.  NUMBER  OF  RAGES 

47 


16  PRICE  COOE 


17.  SECURITY  CLASSIFICATION 
OF  REPORT 

Unclassified 


16.  SECURITY  CLASSIFICATION 
OF  THIS  MGE 
Unclassified 


19.  SECURITY  CLASSIFICATION 

°F  WBStkd 


20.  LIMITATION  OF  ABSTRACT 


MSN  7540-01-280  5500 


i 


Standard  form  298  2  891 

•’•WWO  Or  *r»Vl  MO  .’19  9 


Approved  for  public  release;  distribution  is  unlimited. 


CROSS-SHORE  SEDIMENT  TRANSPORT  ON  A 
NATURALLY  BARRED  BEACH 

by 

Randall  T.  Humiston 
Lieutenant,  United  States  Navy 
B.S.,  University  of  Illinois  at  Carbondale,  1984 

Submitted  in  partial  fulfillment 
of  the  requirements  for  the  degree  of 

MASTER  OF  SCIENCE  IN  PHYSICAL  OCEANOGRAPHY 
AND  METEOROLOGY 

from  the 


Author : 


Approved  by: 


NAVAL  POSTGRADUATE  SCHOOL 
March  1993 


Randall  T.  Humiston 


Edward  B.  Thornton,  Thesis  Advisor 


0. 


Thomas  C.  Lippmann,  Second  Reader 


Curtis  A.  Collins,  Chairman 
Department  of  Oceanography 


Accesion  For 


-2 


NTtS  CRA4I 
OTIC  TAB  0 

Unannounced  0 

Justification 


111 


By 


Distribution/ 


Availability  Codes 


Dist 


i  Avail  and  I  or 


Special 


ABSTRACT 


Bailard's  sediment  transport  model  (1981)  is  evaluated 
using  field  data  obtained  on  a  naturally  barred  beach. 
Principal  field  measurements  consist  of  a  cross-shore  array  of 
bi-directional  current  meters  spanning  the  surf  zone  and  daily 
bathymetric  surveys.  The  model  predicts  bed  and  suspended 
load  transport  separately  based  on  various  velocity  moments. 
The  velocities  are  partitioned  into  mean  currents,  long  waves 
(<  0.05  Hz)  and  short  waves  (>  0.05  Hz)  to  determine  their 
relative  importance  to  the  transport.  Velocity  moments  are 
then  computed  over  90  minute  intervals  to  resolve  tidal 
fluctuations.  Finally,  predicted  transport  rates  are 
integrated  and  compared  with  daily  cross-shore  bathymetric 
profiles  (averaged  over  a  400m  length  of  beach) . 

Results  indicate  that  suspended  load  was  consistently 
greater  than  bed  load,  as  much  as  an  order  of  magnitude  during 
episodes  of  large  incident  waves,  owing  to  the  slow  fall 
velocity  (2cm/s)  of  the  fine  grain  sand  within  the  surf  zone. 
The  contribution  by  the  mean  current,  long  and  short  waves  to 
the  cross-shore  transport  were  of  the  same  order.  Variance  of 
transport  during  all  stages  of  tide  and  over  a  range  of 
incident  waves  were  consistently  greater  in  the  vicinity  of 
the  bar  and  trough  than  seaward  of  the  bar  and  on  the  beach 
face.  Tidal  signatures  were  apparent  in  all  modes  of 
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transport.  The  model  appeared  to  under-predict  measured 
bathymetry  during  low-energy  periods  and  over-predict  during 
high-energy  conditions.  However,  the  model  does  correctly 
predict  the  first  order  movement  of  the  bar. 
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I. 


INTRODUCTION 


The  idea  that  bed  load  transport  rate  of  sediment  grains 
in  a  fluid  can  be  related  to  the  energy  expended  on  transpor¬ 
ting  the  grains  by  the  fluid  flow  in  a  stream  was  introduced 
by  Bagnold  (1956).  He  further  extended  this  hypothesis  to 
include  the  relationship  between  the  interior  flow  and  the 
suspended  load  transport  rate  (Bagnold  1963,  1966).  As  the 
grains  are  at  different  heights  in  the  water  column  for  the 
two  modes  of  transport  the  forces  acting  on  the  grains  are 
distinctly  different.  A  mode  dependent  efficiency  factor  was 
introduced  to  account  for  this  difference. 

Bowen  (1980)  and  Bailard  (1981)  independently  developed  a 
transport  model  for  oceanic  beaches  based  on  the  work  of 
Bagnold.  Bailard' s  model  retained  only  the  first  two  terms  of 
the  expansion  of  the  velocity  moments  and  assumes  no 
autosuspension.  The  model  predicts  the  total  (i.e.,  bed  and 
suspended  load)  longshore  and  cross-shore  transport  rate  over 
an  arbitrary  bottom  topography.  Guza  and  Thornton  (1985) 
analytically  examined  the  relative  importance  of  the  velocity 
moments  in  Bailard 's  model  using  monochromatic  and  Gaussian 
velocity  distributions,  and  compared  these  results  with  field 
measurements.  They  found  modeling  the  nearshore  velocity 
field  as  a  linear,  Gaussian,  random  process  gave  reasonable 
predictions  for  some  moments  and  poor  predictions  of  others. 
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Monochromatic  waves  gave  poor  predictions  of  the  moments. 
They  further  fou  i  that  mean  flow  and  low  frequency  infra¬ 
gravity  waves  must  be  included  to  properly  describe  the 
velocity  field  for  predicting  sediment  transport. 

The  principal  aim  of  this  paper  is  to  evaluate  the  ability 
of  Bailard's  model  to  predict  sediment  transport  from  velocity 
data  measured  on  a  naturally  barred  beach.  The  data  were 
acquired  during  the  nearshore  experiment  Delilah  and  are 
composed  of  eight  closely  space  velocity  measurements 
extending  from  the  shoreline  to  across  the  bar.  The  data  are 
unique  in  that  they  were  continuously  recorded  over  a  period 
of  three  weeks  allowing  examination  of  tidal  effects.  Since 
sediment  transport  was  not  measured  directly,  verification  is 
by  comparing  predicted  bathymetric  changes  to  observed 
changes.  A  secondary  goal  is  to  determine  the  relative 
importance  of  forcing  by  mean  currents,  long  waves  and  short 
waves  under  changing  conditions  over  tidal  cycles. 
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II.  MODEL 


Bailard's  (1981)  time  averaged,  immersed  weight  sediment 
transport  rate,  <1(0  >  ,  is  given  by 


<1(0  >  =  pCf 


tan<J> 


< | u (  t)  !JU(C)>  -  ||£|<|u(t)  !3>  1 


(1) 


<|u(  t)  |3u(  t)  >  -  tanp  <|u(  t)  |5>  1 

It 


where  ju(  t)  {  *  [Q2  +  P2  *  TP  +  v2  *  2  (tiu  +  Vv)  ]  2  is  measured  at  the 
top  of  the  bottom  boundary  layer,  p  is  the  water  density,  Cf 

is  the  bed  drag  coefficient,  p  is  the  beach  slope,  <J>  is  the 
internal  friction  angle  of  the  grains  in  the  bed  where 
tan<|>  =  0 .63  ,  and  W  is  the  fall  velocity.  The  vector  I  is  in 
the  +x  (offshore)  direction.  In  equation  (1)  the  terms  in  the 

first  sot  of  brackets  represent  the  bed  load  contribution  to 

total  sediment  transport  and  the  terms  in  the  second  set  of 

brackets  represent  the  suspended  load  contribution.  The 

suspended  and  bed  load  efficiency  factors,  €,  and  €b  are  the 
ratio  of  stream  power  to  suspended  and  bed  load  work  rates, 
respectfully. 

Following  Guza  and  Thornton  (1985),  the  bi-directional 
velocity  components  are  partitioned  into  oscillatory  (tilde) , 
and  mean  (overbar)  flows 
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fif(t)  =  (0  +  u)  l  -*•  (tf  ♦  v)  j  (2) 

Substitution  of  equation  (2)  into  equation  (1)  gives  equations 
for  time  averaged  alongshore  and  cross-shore  sediment 
transport  rates.  The  equation  of  interest  for  our  work  is  the 
cross-shore  transport 


<I{x)>-pCe 


tan4» 

s 


»  2  3  4 

<Q3>  *■  <QP3>  +  <Q2  +  P3>lj+2<n2>ij 


61  s 

+  2  <£!?>  v+tJv3  +  u1  -  -~anr  <  |d(  t)  |3> 

tan$  ' 


*  10  it 

<|0(t>  |Jfl>  +  |3u>  -_ftanp<|iJ(t)  |5> 

It 


(3) 


The  oscillatory  terms  (0,  V)  are  further  expanded  into 
infragravity.  0j  (f<0.05  Hz),  and  incident,  0,  (f>0.05  Hz), 
wave  components  0(  t)  =  Qj  ( t)  +  Q0(  t)  .  Cross-shore  transport  was 
shown  to  be  insensitive  to  the  sign  of  longshore  current. 
Determining  the  individual  contribution  by  mean  current,  long 
and  short  wave  is  accomplished  by  allowing  the  contribution  of 
interest  and  setting  the  other  terms  to  zero  (e.g.,  to  find 
transport  due  only  to  the  long  wave  component,  Uj  and  vj  are 
as  measured  and  us,  v5,  u,  v  «  0) .  This  procedure  eliminates 
cross  terms  that  normally  contribute  to  the  transport. 

In  addition  to  partitioning  the  contributions  to  the 
sediment  transport  due  to  mean  currents  and  oscillatory  long 
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and  short  waves ,  the  relative  magnitudes  of  the  eleven  terms 
of  equation  (3)  are  also  determined. 

The  volume  sediment  transport  per  unit  width,  q,  (King  and 
Seymour,  1986)  is  calculated  from  the  immersed  weight  sediment 
transport  per  unit  width,  i 


<7 


i 

(PM  -  p)  9 


(4) 


where  p%  -  2.65  gm/cm3  is  the  density  of  quartz  sand.  To 
predict  changes  in  bathymetry  due  to  cross-shore  transport  it 
is  assumed  the  bottom  contours  are  straight  and  parallel  and 
conditions  are  homogeneous  alongshore  such  that  *  0 . 

Conservation  of  mass  is  used  to  obtain  an  average  change  in 
depth  between  current  meters  due  to  cross-shore  sediment  flux 


dh 

Tt 


(5) 


where  n  -  0.7  is  introduced  to  account  for  packing  of  the 
loose  grains  in  the  bed,  and  h  is  depth  in  meters. 
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Ill .  DELILAH  EXPERIMENT 


Data  were  acquired  as  part  of  the  DELILAH  nearshore 
processes  experiment  at  the  Coastal  Engineering  Research 
Center's  Field  Research  Facility  (FRF)  in  Duck,  NC  during 
October  1990.  The  period  investigated  is  from  the  morning  of 
8  October  until  the  morning  of  14  October.  The  environmental 
conditions  are  shown  in  Figure  1.  This  period  was  chosen 
because  it  contains  both  periods  of  relatively  calm  seas  and 
storm  conditions.  Of  note  is  the  increased  wind  speed, 
significant  wave  height  (HMO)  and  alongshore  current  during 
the  storms  commencing  on  10  October  and  ending  of  the  13th. 
On  the  10th,  a  frontal  system  from  the  south  arrived  resulting 
in  waves  up  to  2m  incident  at  relatively  large  angles  (about 
40  degrees  from  the  south)  driving  strong  northward  longshore 
currents  (up  to  1.5  m/s).  On  the  13th,  narrow  band  swell 
waves  up  to  2.5m  arrived  from  hurricane  Lili.  Although  these 
waves  were  larger,  the  incident  wave  angles  were  less  and  the 
resulting  longshore  currents  were  not  as  large. 

The  minigrid  is  defined  as  the  area  of  bathymetry  surveyed 
daily  by  the  FRF's  staff  with  the  Coastal  Research  Amphibious 
Buggy  (CRAB) .  The  bathymetry  before  the  large  incident  waves 
on  11  October  shows  significant  variability  alongshore  and 
across-shore  (Figure  2) .  The  averaged  (over  400m)  bathymetry 
profiles  during  the  experiment  are  shown  in  Figure  3.  A  tidal 
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plateau  existed  until  the  start  of  the  storms  when  a  well 
developed  bar  was  formed.  After  the  storms  commencing  on  11 
October  much  of  the  alongshore  variation  has  vanished  and  a 
linear  longshore  bar  is  present  for  the  remainder  of  the 
experiment. 

The  velocity  data  used  were  obtained  from  a  cross-shore 
array  of  bi-directional  current  meters  located  at  about  the 
center  of  the  DELILAH  minigrid  (Figure  4)  .  Also  shown  in 
Figure  3  are  the  location  of  current  meters.  The  cross-shore 
distance  in  meters  relative  to  the  FRF  coordinate  system  and 
depth  relative  to  National  Geodetic  Vertical  Datum  (NGVD)  are 
given  in  Table  1.  Cross-shore  distance  and  depths  are  in 
meters.  Data  from  the  most  shoreward  of  current  meters  in  the 
array,  cmlO,  were  not  used  as  the  sensor  was  frequently  out  of 
the  water.  At  about  0700  on  the  11th,  cm60  was  lost.  From 
the  11th  on,  bathymetry  predictions  are  averaged  between  cm50 
and  cm70  versus  between  cm50  to  cm60  and  from  cm60  to  cm70. 

Velocity  component  data  were  numerically  rotated  to  have 
axes  parallel  and  perpendicular  to  depth  contours.  An  average 
beach  orientation  was  determined  for  each  day  based  on  the  2m 
contour  averaged  over  about  400m  alongshore.  The  rotation  is 
such  that  the  alongshore  direction  is  defined  in  the  direction 
of  the  2m  contour.  The  rotations  for  each  day  are  given  in 
Table  2.  The  wave  induced  horizontal  velocities  were  depth 
corrected  from  sensor  depth  to  the  top  of  the  bottom  boundary 
layer  to  conform  to  the  requirements  of  Bailard's  model  using 
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a  transfer  function  based  on  linear  wave  theory,  which  is 
generally  a  good  approximation 


H(f)  * 


_ 1 

coshJc(h  ♦  zm) 


(6) 


where  k  is  the  wave  number,  h  is  the  depth  and  is  the 
current  measurement  depth.  Data  from  tide  gage  511,  located 
outside  the  surf  zone  near  the  8m  depth  contour,  were  used  to 
calculate  the  mean  water  level  relative  to  NGVD  to  determine 
the  90  minute  mean  h  in  equation  6. 

Typical  energy  density  spectra  are  shown  in  Figure  5.  It 
can  be  seen  that  there  is  a  spectral  gap  at  about  0.05  Hz 
which  makes  it  a  logical  frequency  to  partition  the  data  into 
long  and  short  waves.  This  was  done  in  order  to  determine  the 
influence  of  different  classes  of  waves  on  sediment  transport. 

Sediment  ranged  in  size  from  about  0.12-0. 2mm  within  the 
surf  zone  and  consisted  mostly  of  nearly  spherical  quartz. 
Cross-shore  size  distribution  is  shown  in  Table  3. 
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IV.  APPLICATION  OP  MODEL 


In  the  application  of  Bailard's  model  to  field  data,  the 
efficiency  factors  €s  and  and  the  fall  velocity,  W,  are 
specified  based  on  mean  grain  size.  The  fall  velocity  for 
spherical  quartz  grains  of  mean  size  0.2mm  is  2  cm/sec  (Dyer, 
1986).  A  nominal  €a  =  0.015  for  a  stream  is  used  as  determined 
by  Bagnold  (1966)  and  is  used  by  others  (e.g.,  Bowen  1980, 
Bailard  1981) .  Bagnold  (1966)  presents  a  functional 
relationship  for  €b  with  stream  velocity  and  grain  size. 
Using  the  mean  grain  size  range  of  0.12-0. 2mm  found  at  Duck 
(Table  3),  and  stream  velocities  with  range  of  1- 
2m/sec,  €6»0.135  ±0.004  ,  which  indicates  eb  is  not  sensitive 
to  the  observed  range  of  grain  sizes  or  velocities. 

A  value  of  Cf  *  0.003  was  used  which  is  the  value  obtained 
from  longshore  current  modeling  of  the  DELILAH  data  (Church 
and  Thornton,  1993) . 

The  values  for  tan£  were  calculated  from  the  measured 
bottom  slopes  in  the  vicinity  of  each  current  meter. 

The  total  immersed  weight  sediment  transport  partitioned 
into  bed  and  suspended  load  components  for  current  meters 
cm20-50  and  cm80  are  shown  in  Figures  6-10.  Total  cross-shore 
transport  is  simply  the  sum  of  bed  and  suspended  load.  Tidal 
variation  in  depth  is  superimposed  on  all  graphs  for 
comparison  (tidal  amplitudes  are  not  to  scale) .  Negative 
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transport  indicates  periods  of  onshore  transport  (accretion) 
and  positive  transport  indicates  offshore  transport  (erosion) . 
Times  of  high  energy  (Figure  1)  are  easily  discernable  at  each 
current  meter  location.  Fluctuations  in  transport  are  roughly 
of  the  same  period  as  the  tide,  although  magnitude  and  phase 
vary  with  time  and  location.  Transport  due  to  suspended  load 
is  consistently  greater  than  bed  load,  particularity  during 
times  of  high  waves  (11  October,  and  13  October) . 

Also  shown  is  the  transport  due  to  mean,  long  wave  and 
short  wave  components  of  the  flow,  also  partitioned  into  bed 
and  suspended  load.  As  can  be  seen  from  equation  (3)  ,  the 
components  of  the  flow  contribute  non-linearly  to  the 
transport,  and  thus  cannot  merely  be  summed  to  equal  the 
total. 

The  dominant  forcing  varies  with  location  and  conditions 
and  will  be  discussed  in  order  of  increasing  distance  from  the 
beach.  In  the  area  between  the  foreshore  and  trough  (near 
cm20)  there  is  very  little  accretion,  with  most  of  the  erosion 
occurring  during  storm  conditions.  In  this  region,  cross¬ 
shore  transport  is  dominated  by  longshore  and  cross-shore  mean 
flow.  The  short  waves  have  a  slightly  greater  influence  than 
the  long  waves. 

In  the  trough  (Figure  7,  cm30)  the  major  events  are 
accretionary.  Short  wave  flow  dominated  the  total  transport 
in  the  trough  during  the  second  storm  on  13  October. 
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On  the  shoreward  slope  of  the  bar  (Figure  8,  cm40)  both 
accretion  and  erosion  were  evident  with  short  wave  component 
accretxon  usually  dominating.  Of  note  is  that  the  long  wave 
component  produced  mostly  onshore  transport. 

On  top  of  the  bar  (Figure  9,  cm50)  transport  was  almost 
exclusively  offshore  with  mean  and  short  wave  flow  nearly 
equal  and  both  greater  than  long  wave  flow.  In  this  area  the 
periods  of  high  transport  rates  were  of  greater  duration  than 
in  any  other  area.  A  definite  tidal  period  is  seen  here  with 
mean  and  long  wave  flow  contributions  being  predominantly 
during  low  tide  and  the  short  wave  contribution  being  during 
transition  from  both  low  to  high  and  high  to  low  tide  when  the 
tidal  currents  can  be  expected  to  be  at  their  maximum. 

Further  seaward  at  cm80  (Figure  10)  transport  is  dominated 
by  the  short  waves,  and  is  confined  to  periods  of  very  high 
waves . 

Of  particular  interest  are  two  cases  in  which  the 
suspended  and  bed  load  transports  are  180*  out  of  phase.  The 
first  is  at  cm20  for  the  short  wave  flow  component  on  the 
evening  of  the  13th.  The  second  is  at  cm40  for  the  mean  flow 
component  at  roughly  the  same  time.  The  total  suspended  and 
bed  load  transport  rates  for  these  cases  are  in  phase. 

The  contributions  to  the  total  transport  by  each  of  the 
eleven  terms  of  equation  (3)  were  evaluated  by  the  model  as 
well.  As  suspended  load  was  calculated  to  be  much  greater 
than  bed  load,  only  one  representative  plot  of  the  three 


11 


r 


suspended  load  terms  is  shown  (Figure  11) .  The  dominance 
seems  to  vary  between  terms  9  and  11,  depending  on  time  and 
location.  All  three  terms  exhibited  periods  of  positive  and 
negative  values  and  differing  phase  relationships  between  each 
other. 

The  daily  predicted  change  in  depth  between  current  meter 
positions  was  determined  using  equations  (4)  and  (5)  (Figure 
12)  and  compared  with  the  observed  bathymetric  profiles.  In 
the  application  of  equation  (5)  to  calculate  the  changes  in 
the  bathymetry,  the  sediment  fluxes,  calculated  every  90 
minutes  are  integrated  (summed)  between  times  of  the 
bathymetry  surveys  (approximately  every  24  hours)  such  that 

Ah  "  c)  c>]  (7) 

where  At  is  the  time  between  profile  transects  (nominally 
24hrs)  and  Ax  =  Xj+1-  the  distance  between  current  meters. 
An  average  profile  was  used  to  minimize  the  effect  of 
longshore  inhomogeneity.  For  each  day,  bathymetry  profiles 
were  averaged  over  200m  north  and  200m  south  of  the  current 
meter  array. 

There  is  little  discernable  difference  between  predicted 
and  observed  depths  from  8-9  October.  A  trough  begins  to  form 
between  cm20  and  cm30  from  9-10  October.  The  model  does  not 
accurately  account  for  this  change.  Between  the  10th  to  llth 
the  trough  continues  to  erode  at  an  accelerated  rate  and  the 
bar  migrates  to  seaward.  The  model  predicts  a  decrease  in 
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depth  between  cm20  and  cm30,  opposite  of  the  observed 
hollowing  of  the  trough.  The  model  does,  however,  indicate  a 
deepening  of  the  trough  between  cm30  and  cm50.  Between  cm50 
and  cm60  the  model  correctly  predicts  the  migration  of  the 
bar.  During  the  low  waves  between  the  11th  and  12th  there  is 
little  change  in  bathymetry,  also  reflected  by  the  model.  As 
the  waves  increase  during  the  12th  to  13th  the  trough  deepens 
and  the  bar  migrates  further  seaward.  The  model  again 
incorrectly  shows  a  shallowing  between  cm20  and  cm30,  but 
correctly  shows  the  trough  deepening  and  seaward  bar 
migration.  From  the  13th  to  14th  the  model  predicts  much  more 
transport  than  actually  occurs.  Interestingly,  the  transport 
is  consistent  with  those  days  which  did  have  significant 
change,  (i.e.,  deepening  of  the  trough  and  seaward  migration 
of  the  bar) . 
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V.  DISCUSSION 


Suspended  load  dominates  bed  load  owing  to  the  slow  fall 
velocity  (0.02m/s)  of  the  fine  grain  sand  (Table  1).  This 
slow  fall  velocity  means  that  the  grains  remain  in  suspension 
much  longer  with  less  turbulent  energy  required. 

To  a  first  order  approximation,  the  model  correctly 
predicts  the  deepening  of  the  trough  and  offshore  migration  of 
the  bar.  However,  in  general,  it  under  predicts  changes  in 
depth  during  low  energy  wave  periods  and  over  predicts  during 
high  energy  episodes. 

Between  the  foreshore  and  trough,  transport  is  dominated 
by  mean  flow.  The  strong  longshore  current  exerts  a  strong 
bottom  stress  while  the  weaker  cross-shore  flow  (undertow) 
determines  the  direction  of  transport.  The  model  consistently 
predicts  a  convergence  in  transport  between  cm20  and  cm30, 
resulting  in  a  decrease  in  depth.  This  is  contrary  to 
observed  bathymetry  changes.  One  possibility  is  that  sediment 
is  being  transported  out  of  the  system  by  the  longshore 
current.  While  only  calculating  the  cross-shore  component  of 
sediment  transport,  one  must  assume  that  the  longshore 
transport  is  constant  and  that  any  longshore  variation  in 
morphology  can  be  averaged  to  a  spatially  constant  mean 
profile.  Averaging  adjacent  profiles  may  not  be  sufficient  to 
account  for  divergences  and  convergence  in  the  longshore 
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drift.  As  pointed  out  earlier,  the  bathymetry  at  the 
beginning  of  the  experiment  was  significantly  3D  (longshore 
and  cross-shore  variations)  but  evolved  into  a  2D 
configuration  with  a  significant  longshore  bar  after  the  storm 
on  the  11th,  Figure  2.  This  suggests  that  the  assumption  of 
longshore  homogeneity  was  not  always  valid.  Using  only  the 
cross-shore  component  of  transport,  particularly  during  times 
of  3D  bathymetry,  may  not  be  sufficient  to  account  for  the 
observed  changes  in  the  bathymetry. 

The  calculated  divergence  of  the  sediment  flux  was  not 
well  resolved  with  the  current  meters  separated  by  approxi¬ 
mately  20m.  Inspection  of  Figure  3  shows  that  cm20-50 
positioning  meant  that  the  slope  of  the  bottom  changes  sign  at 
nearly  each  of  these  sensors.  As  term  11  of  equation  (3)  is 
a  function  of  this  slope  and  one  of  the  dominant  terms,  the 
calculated  sediment  flux  for  term  11  likewise  changed  sign  at 
each  current  meter. 

The  model  assumes  that  current  velocities  are  taken  near 
the  top  of  the  bottom  boundary  layer.  The  wave  induced 
velocities  were  depth  corrected  using  the  linear  wave  theory 
transfer  function  (equation  (6)).  However,  the  mean  currents 
were  not  depth  corrected  and  substantial  error  could  be 
introduced  when  considering  the  undertow. 

Transition  time  from  quiescence  to  high  energy  episodes  as 
seen  in  Figure  1  is  relatively  short.  During  quiescent  times 
the  main  influence  on  transport  and  bathymetry  changes  is  by 
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the  tide.  Daily  bathymetric  transects  are  aliased  over  these 
times.  Reducing  the  time  between  surveys  to  several  times 
daily  would  provide  bathymetry  better  suited  to  evaluate  the 
model . 

It  was  not  possible  to  attribute  the  relative  dominance 
between  terms  9,  10,  and  11  of  equation  (3)  to  any  location  or 
environmental  condition  as  the  terms  varied  greatly  in  time 
and  space  as  well  as  relative  phase.  The  only  consistent 
feature  was  a  definite  tidal  signature.  That  the  tide  is 
apparent  in  the  model  transport  rates  indicates  that  future 
measurements  of  sediment  flux  should  be  taken  continuously 
over  tidal  cycles. 
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VI .  CONCLUSIONS 


Results  indicate  that  suspended  load  was  consistently 
greater  than  bed  load  during  the  DELILAH  experiment,  as  much 
as  an  order  of  magnitude  during  episodes  of  large  incident 
waves,  due  to  the  small  grain  size  and  resultant  slow  settling 
velocity. 

The  three  velocity  moments  for  suspended  transport  varied 
both  spatially  and  temporally  making  evaluation  of  the 
dominant  forcing  components  inconclusive. 

The  mean,  long,  and  short  wave  contribution  to  the  cross¬ 
shore  sediment  transport  were  all  of  the  same  order.  The  mean 
current  forcing  was  dominated  by  the  strong  longshore  current 
during  the  experiment  with  the  cross-shore  current  at  the  bed 
determining  the  direction  of  the  transport. 

Variance  of  transport  at  all  tide  stages  and  during 
different  wave  conditions  was  much  greater  in  the  vicinity  of 
the  bar  and  trough.  A  tidal  signature  was  evident  throughout 
the  cross-shore  for  each  mode  of  flow  suggesting  that  future 
experiments  should  acquire  sediment  flux  data  continuously 
over  the  tidal  cycle,  or  that  in  predicting  bathymetric 
changes  the  integration  time  step  should  resolve  the  tidal 
cycle. 

Model  predictions  of  bathymetric  changes  appear  to  under¬ 
predict  measured  bathymetry  during  low-energy  periods  and 
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over-predict  during  high-energy  conditions.  However,  the 
model  does  correctly  predict  the  first  order  movement  of  the 
bar. 

Closer  cross-shore  placement  of  flow  sensors  in  the 
vicinity  of  the  trough  and  bar  and  several  bathymetry 
measurements  per  day  could  result  in  more  accurate 
representation  of  the  transport  processes  and  evaluation  of 
the  model. 
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APPENDIX  A 


TABLE  1.  CURRENT  METER  cnl0-90  CR088-8H0RE  DISTANCE 

IN  FR7  COORDINATES  AND  DEPTH  RELATIVE  TO  NGVD 


Current 

Meter 

Cross-shore 
distance  (m) 
FRF 

Depth  (m) 
NGVD 

cmlO 

125 

o 

• 

o 

-4 

cm20 

145 

-0.77  | 

cm30 

170 

-.066  | 

cm40 

189 

-.082 

cm50 

207 

-.098 

cm60 

226 

-1.01 

cm70 

245 

-1.05 

cm80 

295 

.  “1*66 

cm90 

370 

-1.68  1 

TABLE  2.  ROTATION  ANGLE  FOR  DEFINING 
AL0NG8H0RE  DIRECTION 


Date 

Rotation  angle 

Oct  8 

3.0 

Oct  9 

4.4 

Oct  10 

4.0 

Oct  11 

2.5 

Oct  12 

2.6 

Oct  13 

2.5 

Oct  14 

1.1 
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TABLE  3.  SEDIMENT  GRAIN  SIZE  ACROSS- SHORE 


Cross-shore 
distance  (ra)  FRF 

Sediment  mean 
diameter  (mm) 

f  138 

.122 

f  149 

.182 

3  192 

.205 

256 

.143 

328 

.142 

484 

.158 
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Figure  2.  Minigrid  3D  plot  of  bathymetry  on  9  October  (top)  and 
11  October  (bottom)  (from  Birkemeir,  Hathaway,  Smith,  Baron, 
Leffler,  and  many  others,  1991). 
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Figure  4.  DELILAH  minigrid.  Longshore  and  cross-shore  distances 
are  in  FRF  coordinates  (from  Berkemeir,  Hathaway,  Smith,  Baron, 
Leffler,  and  many  others,  1991) . 
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Figure  6.  Immersed  weight  sediment  transport  per  unit  width  at  cm20.  Upper  left  is  the 
total,  upper  right  due  to  mean  currents,  lower  left  due  to  long  wave  component  of  curren 
lower  right  due  to  short  wave  component  of  current.  Each  transport  is  partitioned  into 
suspension  and  bed  load.  Tides  are  superimposed  and  amplitude  is  not  to  scale. 
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7.  Immersed  weight  sediment  transport  per  unit  width  at  cra30.  Upper  left  is  the 
upper  right  due  to  mean  currents,  lower  left  due  to  long  wave  component  of  curren 
right  due  to  short  wave  component  of  current.  Each  transport  is  partitioned  into 
sion  and  bed  load.  Tides  are  superimposed  and  amplitude  is  not  to  scale. 
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Figure  10.  Immersed  weight  sediment  transport  per  unit  width  at  cm80.  Upper  left  is  the 
total,  upper  right  due  to  mean  currents,  lower  left  due  to  long  wave  component  of  current 
lower  right  due  to  short  wave  component  of  current.  Each  transport  is  partitioned  into 
suspension  and  bed  load.  Tides  are  superimposed  and  amplitude  is  not  to  scale. 
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Figure  11.  Model  terms  9  thru  11  for  cm50.  Term  9  due  to  to  ociliatory  flow,  term  10 

due  to  mean  flow  and  term  11  due  to  bottom  slope.  Tides  are  superimposed  and 
amplitude  is  not  to  scale. 
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Figure  12.  Predicted  changes  in  bathymetry.  Circles  are 
predicted  change  over  24  hours  to  initial  observed  profile,  solid 
line.  Dotted  line  is  observed  profile  24  hours  later.  Given  a 
perfect  prediction,  the  circles  would  fall  on  top  of  the  dotted 
li 
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